LASTiC SCATTE 
!LE PION PROD1 
AT 0.939 Bev y 








DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



UCRL-8302 

Physics and Mathematics 

UNIVERSITY OF CALIFORNIA 

Radiation Laboratory 
Berkeley; California 

Contract No. W-7405 -eng - 48 



ir“-p ELASTIC SCATTERING AND SINGLE PION PRODUCTION 

AT 0.939 Bev/c 

Lee Baggett, Jr. 

May 28; 1958 



35 /• 



Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE 
IN 

PHYSICS 

United States Naval Postgraduate School 
Monterey, California 



Printed for the U S. Atomic Energy Commission 



Mir’- AvXG vf. 















































tt"-p elastic scattering and single pion production 

AT 0.939 Bev/c 

by 

Lee Baggett, Jr. 



This work is accepted as fulfilling 
the thesis requirements for the degree of 

MASTER OF SCIENCE 
IN 

PHYSICS 
from the 

United States Naval Postgraduate School 



it - p ELASTIC SCATTERING AND SINGLE PION PRODUCTION 

AT 0.939 Bev/c 

Contents 

Abstract ........ iv 

I. Introduction . 1 

II. Experimental Procedure . 3 

III. Scanning Efficiency and Acceptance Criteria 

A, Scanning Efficiency 

General 3 

Depth Dependence 3 

x-y Dependence 4 

Orientation of Outgoing Prongs 8 

Repeated Scans . 8 

B. Acceptance Criteria 8 

IV. Track Analysis 

A. Spatial Reconstruction of a Track 10 

B. Measurement Errors . 13 

Internal Measurement Errors .13 

Turbulence Errors 15 

Incident - Beam Momentum ZZ 

V. Kinematic Analysis 

A. General . . Z4 

B Computation of Measurement Errors and Coplananty . Z4 

C. Graphical Analysis of Inelastic Events 

Ellipse Plots Z6 

Transformation to Center -of -Mass System 31 

D. Computer Program for Kinematic Analysis of Two Prong 

Events 

General 36 

Choice of Dynamical Variables 36 

Interpretation Defined 37 



li 



Kinematic Equations of Constraint . 38 

Definition of a Kinematic Fit 39 

Statistical Estimates of the Quality of the Fit . . . .40 

Trajectory of the Missing Neutral Particle .... 41 

Mathematics of the Fitting Procedure 41 

Statistical Correlation of Adjusted and 

Measurable Variables 43 

Comparison of Measured and Computed Values . 44 

VI. Classification of Events . 49 

VII. Results 

A. Cross Sections .52 

B. Elastic Differential Cross Section 52 

C. Distribution of Charged Prongs in SingLe -Pion 

Production 55 

References 59 



lii 



. 

. 

•' '• 

Iff .1 0*)|M £' IV 

Hi-Wli, n -T • I .Hi d- 3 



tt'-p elastic scattering and single pion production 

AT 0.939 Bev/c 

Lee Baggett, Jr. 

Radiation Laboratory 
University of California 
Berkeley, California 

May 28, 1958 
ABSTRACT 

An investigation of 7r~-p elastic scattering and singie-pion pro- 
duction at 0.939 Bev/c is reported. Five-hundred and forty two-prong 
interactions in the 10 -inch UCRL liquid -hydrogen bubble chamber have 
been classified. Emphasis is given the methods used in the kinematic 
classification of events. The total cross section for tt -p interactions 
is found to be 53.3 ± 2.4 mb. The elastic-scattering cross section is 
21,0 ±1.1 mb. The elastic differential cross section displays a 
narrow diffraction peak in the forward direction. In the center-of- 
mass system, 35.6 % of the elastic events have pion scatters into the 
backward hemisphere. 

+ 0 

In single -pion production the ratio of the number of tt to tt 
productions is 2.49. The momentum distributions of the charged 
particles do not correspond to the phase-space distributions predicted. 
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tt~ - p ELASTIC SCATTERING AND SINGLE PION PRODUCTION 

AT 0.939 Bev/c 

Lee Baggett, Jr. 

Radiation Laboratory 
University of California 
Berke ley, California 

May 28, 1958 

I. INTRODUCTION 

An experimental study of elastic scattering and single -pion - 
production interactions of negative pions of 0.939 Bev/c momentum in 
hydrogen is reported here. 

Pion proton interactions in this Bev region leading to associated 

1 Z 3 

production have been investigated recently. ; * Inelastic pion reactions 

not involving "strange particles" have been neglected except for early, 

4 5 

low-statistics experiments using diffusion cloud chambers , and 
nuclear emulsions. ^ Elastic pion^proton scattering in this energy 

region, on the other hand, has received additional treatment using a 

7 \ 

propane chamber and, more recently, liquid hydrogen bubble chambers. 

The study described here is based on data obtained as a by- 
product of a recent as sociated -production experiment with the 10-mch 

hydrogen bubble chamber at the University of California Radiation 

9 

Laboratory. Approximately 600 two-prong pion events have been 
measured at each of three (nominal) beam momenta: 0.937, 1.12, and 

1.24 Bev/c. The preliminary kinematic analysis of the 0.927-Bev/c 
data is reported here. Analysis of the 1.12° and 1.24-Bev/c momenta 
data will be reported later. 

The following three reactions are considered: 

TT + p TT + p, 

- , , 0 

TT +p-*7T t p T TT , 

7 T +p-*7T +TT + + n. 

The decision to limit the kinematic analysis to the three listed reactions 
assumes that the frequency of events with the production of two or more 
secondary pions is negligible. The rarity of four -prong events, 
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tt 4-p-^TT +p + Tr + tt, supports this assumption. In principle, 
measurement of the charged -particle tracks in a magnetic field, and 
conservation of momentum and energy, determine unambiguously 
whether one or more neutral particles are involved. In practice the 
identification of missing neutral particles is made more difficult by 
measurement inaccuracies. Classification of multiple -pion-pr oduction 
events requires well-established criteria as to what constitutes a good 
kinematic fit in single -pion production. An approximate upper bound 
for multiple pion production will be established later. 

This report treats the interpretation of the experimental 
results in only a cursory way. Specification of a single -pion-pr o~ 
duction event requires five intrinsic parameters, in contrast with the 
single center -of -mas s scattering angle required for an elastic event. 
Consequently a detailed interpretive treatment of the data is relegated 
to the future. To keep the experiment in perspective, however, it is 
useful to list objectives of a study of elastic and single -pion-pr oduction 
events in this energy region. 

a. Elastic tt -p scattering 

1. Determination of the size and opacity of the proton. 

2. Evaluation of the relative weights that the various angular- 
momentum states contribute to the interaction. 

3. Test of the validity of the dispersion relations (15% effect). 

4. Confirmation of the previously asserted backward -s cattering 
peak (at higher energies). 

b. Inelastic single -secondary-pion production 

1. Isotopic spin dependence of the interaction; the ratio of 

to n production (a priori of equal phase space). 

3 10 

2. Evidence for the existence of a T = J = nucleon isobar. 

3. Correlation of the inelastic pion production with associated 
production. 
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II. EXPERIMENTAL PROCEDURE 
The it -p interactions analyzed in this experiment were selected 
from a group of photographs taken in the 1 0 -inch-diameter hydrogen 
bubble chamber. The chamber was provided with a 10, 900-gauss 
magnetic field. Details of the hydrogen bubble chamber used have 
been described. ^ The exposure was the cooperative effort of the 
Alvarez group, and was primarily for the purpose of investigating 
associated production of strange particles at 0.927 Bev/c nominal 

12 

beam momentum. The it -beam trajectory is discussed elsewhere. 
About 24,000 Bevatron pulses were photographed during the run. 
Approximately 4,800 of these were used in this analysis. Ail film 
was scanned by the professional scanning staff. 



III. SCANNING EFFICIENCY AND ACCEPTANCE CRITERIA 

A. Scanning Efficiency 



1. General 

At each transit of the tt beam through the bubble chamber, two 
stereophotographs of the chamber were taken. After processing, the 
35 -mm film was scanned for all events of interest. Pion interactions 
were recorded, according to frame number and type, by means of 
check marks in appropriate columns of the scan sheets. Events 
entered on these scan sheets were subsequently examined in detail, 
sketched, and measured. The 13 rolls of film used in this experiment 
were scanned twice. Several factors which reduce the probability of 
seeing an event are discussed, together with their compensating 
corrections . 

2. Depth Dependence 

An event with its vertex near the top or bottom glass plates of 
the chamber is likely to be missed because one or both of the emergent 
prongs can be very short. An analysis of this chamber depth dependence 
was based on a comparison of the vertical distribution of 575 event 
vertices and the vertical distribution of 1126 beam tracks. The chamber 
was considered as divided into five horizontal regions. The median 
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plane of the central region, C, coincides with the median plane of the 
bubble chamber; depth of region C is 3.2 in. Regions A and E are the 
extreme regions, extending to the top and bottom glasses, respectively. 
Regions B and D lie between the central region and regions A and E, 
respectively. These divisions are each 0.8 in. in depth. 

The depths of beam tracks in every 50th frame were measured. 
The vertical position of each event vertex was obtained by measure- 
ment on the precision measuring engine. The ratio of event density 
to track density within each horizontal section of the chamber is a 
measure of scanning efficiency within that region. These ratios are 
as follows. 

Table I 



Chamber Depth Dependence 



Region 


Events (N) 


Beam Tracks(M) 


Ratio (N/M) 


A 


2 


2 


i oo +4 - 80 

1,UU -0.83 


B 


3 


12 


0 25 +0 - 3 ° 
_ Q 1? 


C 


563 


1082 


0 52 +0 ° 4 
-0.04 


D 


6 


22 


0 27 +0 ' 22 
-0.14 


E 


1 


8 


0 12 +0 - 27 
d - 0.12 


Total 


575 


1126 


0.51 



Obviously scanning efficiency for events in depth region E is poor. The 
ratios for ail other regions agree within statistics, and only the events 
found in these four regions were considered. No correction to scanning 
efficiency as a function of chamber depth was made. 

3 . x-y Dependence 

A reaction that occurs near the chamber wail may also be 
difficult to detect because of short tracks. Aside from this consideration, 
momentum and angle measurements become less accurate as the length 
of track available for measurement decreases. Hence, long tracks 
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were desired. A central area was superimposed on one stereo view at 
the stereo projector* using chamber fiducials for positioning; only 
those events whose vertices projected into the area were considered 
for measurement. This criterion not only increased the scanning 
efficiency* and the probability of having long tracks* but it also made 
identification of events from the decay* interaction* or stopping of one 
of the emergent prongs more likely. 

In order to detect any variation in scanning efficiency within 
the acceptance area* this area was divided into five equal regions* as 
shown in Fig. 1. The distribution of 600 projected vertex positions 
was compared with the distribution of 8562 beam tracks crossing the 
area. All tracks in every tent-h frame of the experiment were counted. 
On the histogram of beam track distribution, Fig. 2, the number of 
events within each subdivision of the acceptance area has been indicated. 
The ratio of event density to track density within each region of the 
acceptance area is a measure of the scanning efficiency for events 
within that region. These ratios are as follows. 



Table II 



x~y Dependence 


Region 


Events (N) 


Beam Tracks(M) 


Ratio(N/M) 


i 


91 


1415 


0.064 ± .009 


2 


122 


1784 


0.068 ± .008 


3 


152 


2115 


0.072 ± .007 


4 


141 


1871 


0.075 ± .009 


5 


94 


1377 


0.068 ± .009 


T otal 


600 


8562 


0.070 



These ratios indicate that scanning efficiency is uniform, within 
statistical accuracy, throughout the acceptance area, and that no 
dependence of scanning efficiency on position in this area need be 
considered. 
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Fig. 1. Projected view of bubble chamber, showing acceptance 
area. 
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Fig. 
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2. Projected x-y distributions within acceptance region. 
Histogram of distribution of 8562 beam tracks. (600 events) 
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4. Orientation of Outgoing Prongs 

Events that have one or both prongs near the vertical are more 
likely to be overlooked than events with prongs more nearly horizontal. 
To study the dependence of scanning efficiency on orientation of the 
plane defined by the two prongs, a plot of the number of events as a 
function of the angle § was constructed. Here $ is the angle between 
the horizontal plane and the plane containing the two prongs. This 
plot is shown in Fig. 3. The scattering should be isotropic in <£, and 
the anisotropy observed is a measure of the dependence of scanning 
efficiency on event orientation. A correction factor of 1.18 ± .03 is 
derived from the data illustrated in Fig. 3. The total number of 
events seen must be multiplied by this correction factor to account 
for the scanning -efficiency dependence on 

5. Repeated Scans 

All film used in this experiment was scanned twice by two in- 
dependent observers. Within the acceptance volume there were 922 
events seen by both observers; 33 events were seen by the first scanner 
but not seen by the second; 21 events were seen by the second scanner 
but not by the first. The indicated efficiency of each scanner is 
97% ± 1%. The total number of observable two-prong events, based 
on these scans, is 983 ± 1. The total number of two-prong events, 
after correction for scanning -efficiency dependence on c£, is 1160 ± 31. 

B . Acceptance Criteria 

Each two-prong event within the acceptance area was examined 
at a scanning table before being measured. Acceptance criteria, 
specified to reduce analysis time and measurement errors, and to 
fully utilize the automatic tracking capability of the measuring engine, 
were imposed on each event. These criteria were: (1) The vertex 

was distinct, i. e. , no track not associated with the event fell within 
one track width of the vertex in either projected view. (2) The 
incident track direction was within 3° of the mean incident beam 
direction. (3) The incident track curvature was approximately the 
same as the mean beam curvature. (4) Not more than thirty beam 
tracks were visible in either view. 
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Fig. 3. Orientation of plane of emergent prongs. Folded 
distribution of 550 events. $ = 0° = flat in chamber. 
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IV. TRACK ANALYSIS 
A. Spatial Reconstruction of a Track 
Each of the two views of an event accepted for measurement 
was projected on the screen of the precision measuring engine 
(Franckenstein) . This servo auto -tracking projector takes a series 
of points along each prong of the event in each of the two stereo views. 
The solution of the optical problem and the determination of a set of 
points in chamber space from the Franckenstein measurements have 
been previously described. 

When the space coordinates x^, y^, (i = 1,2, • • • • , n) for each 
of the n measured points along a track have been determined, a least- 
squares fit of a parabola is made to the projection of these points on 
the x-y plane. (Fig. 4. ) The least -square s fit is made on the x/' and 
y. M points following a translation and rotation to the x u , y M coordinates 
used to simplify the scaling problem for the digital computer. The 
x n , y M coordinates are fitted to the parabola 



y. " = a + b x " + 
7 i l 



c x 

l 



m2 



(The event reconstruction described is programmed for the IBM-650 
computer. ) 

Although the projection of the track on the x-y plane is a 

circular arc, a parabolic fit is made rather than a fit to a higher - 

order curve, because at the momenta measured the increase in 

accuracy is not sufficient to justify more elaborate programming. 

When the parabola has been determined, the r oot -mean -square 

deviation cr of the measured points from the fitted curve is ex- 
xy 

pressed by / \ i 

6 . 2 \ 2 



xy \ n-3 

where n is the number of points measured, and 6 is the deviation 

of each projected point from the fitted curve. For n <3, is not 

meaningful, since a parabola can be fitted exactly to three points, and 

therefore no calculation of O is made. 

xy 

The azimuthal angle, between the tangent to the x-y projection 

LS 

of the track and the x axis, is given by $ - + 6 . The value 



LS 

of 6 is obtained from the coefficients 



R 

b* and c 1 and the half-chord 
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Fig. 4. Space coordinate system and track projections into 
planes . 
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length R : 



R = x M - Xj". 

For each track there are two azimuthal angles, <f> i and <j>^, correspond- 
ing respectively to the beginning and end of each track. The value <j>^ 
is used only with the incident beam track, for on this track alone is 
measurement directed toward the event vertex. 

To specify the length of the track, f , and the dip angle X, the 
coordinates are fitted to a straight line in the x M -z M plane. Although 
it is true that each track is a helix, and hence its projection on the 
x H -z n plane is sinusoidal, the momenta were sufficiently high to allow 
a least-squares fit to be made to the straight line z. ,? = d + ex.". 

The root-mean-square deviation, defined by 

f $ 6. 2 \1 



°Z 



n-2 



is then obtained from the coefficients of the straight-line expression. 

As above, 6^ and n are the deviations of the measured points (pro- 
jected into the x ,, z 11 plane) from the fitted curve and the number of 
points measured, respectively. 

With the z dependence of the track now specified, the dip angle 
X and the track length l are given by the expressions 

tan X = - -g- ; 1 = 2(e 2 + R 2 )* . 

Fitted end points are obtained from the least-squares solution. 

The final quantity of interest- -the projected curvature of the 
track K--is determined from the fit in the x-y plane, according to 
K = l/.3Hp, where p = R /2c. The magnetic field strength, H, is 
evaluated at the midpoint of the track. The momentum can be calculated 
from the curvature K and dip angle X by the expression p = (K cos X) 
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B. Measurement Errors 



The automatic tracking feature of the precision measurement 
projector makes positioning of the stage at each point operator -inde - 
pendent. This feature makes it feasible to accurately determine the 
errors associated with measurements of angle and curvature. 

Semi -empirical formulae for the standard deviation errors in 
the curvature K of a track, the azimuthal angle 4>, and the tangent of the 
dip angle, tan X, were constructed from statistical studies of meas- 
urements made both on beam tracks and on tracks of actual events. 

Errors fall into two categories: (a) internal-measurement 

errors, which reflect the accuracy of the measuring device and the 
intrinsic image quality of track; (b) turbulence errors, resulting from 
distortion of tracks by the turbulent liquid hydrogen during the time 
interval between occurrence of the event and the light flash. 

1 . Internal Measurement Errors 

During measurement, track points are taken at approximately 
1 -centimeter intervals along the projection of the track in the x-y 
plane. A fit to these measured points (x., y.; j = 0, 1, ••• n), assuming 

equal standard deviations for the y setting for all points, is specified 

2 

as follows: the x quantity, 



n 



- ^ i 

J=0 



a 



m \ 

i ..*0 

1 i=0 1 J / 



is required to assume a minimum value for the ,f best M trial value of 

the parameters cu(i = 0, 1, •••m). The intrinsic accuracy of the 

once determined, leads directly to the standard deviations °k'- V’ 

and oj x attributive to measurement technique, 
tan X n 

Let us consider the error in the dip tangent. If the z dependence 
of a track is stated as z = Qg + Qj x, with x the distance along the 
projected chord length of the track, the dip tangent is given by 
tanX = dz/dx = The standard deviation is then 



4101* I it** fl» c««>| •* • T» •• n .... 

«Af? •• *r •«* - 

r* M| .«♦*«* j *r 
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or 



O' = a 
tan X a ^ 



O' 



tan X 



4s o z 
R 4 n 




(iv -1) 



where R is one -half the projected chord Length, N is the number of 

points measured, and C^. is a correction factor for dependence on the 

13 

number of points taken. 

The derivation of error in azimuthal angle follows a similar 
pattern. When a fit has been obtained in the x-y plane of the chamber, 
following coordinate rotation, of the form 



y = a.Q + a^x + a^x , 
then examination of Fig. 4 shows: 



tan(<j>.-<j> R ) = dy/dxj x=() = a { . 



For a 10-cm track with momentum of 150 Mev/c, 

this angle decreases rapidly with increasing momentum. 

1 3 

cos(<|>. -4>£ ) = 1 allows the error in <(>. to be written 




WT a c 4> 
rnTn x * N ' 



= 1 4„ 7°, and 
Setting 



(IV -2) 



On the basis of the equation for the fitted curve, as given above, 
the curvature of a track can be defined as 



K = (d 2 y/dx 2 ) [1 + (dy/dx) 2 ] 3 ^ 2 . 

Evaluation of this expression at x = R gives K = 2 cu , where, by 

2 ^ 

choice of coordinate axes, the (dy/dx) contribution can be neglected 
Propagating the error of fit in a, one can approximate the curvature 

13 Vv 

error by 



a 



K 



ML o c K 

r 2 nTn x y N ’ 



(IV -3) 



Both (J and O „ are calculated for the fit achieved. From 
xy Z. 

repeated measurements of individual events and measurements of 
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beam tracks, mean values of 0 0041 cm and 0.0558 cm, respectively, 
were found. Use of these OF and 0 ^ with Eqs„ (IV- 1) to (IV-3), 

B = 305 cm/Bev/c where B = l/(pK), based on the mean 10.9-kilogauss 
magnetic field, results in the following approximate expressions for 
internal measurement errors: 

o* = — [1 +2]sin«t>|] _1 [1 - 0.76/Nl. (IV4) 

tan K rVn 



l 




1.54 

r\Tn 



[1 - l/N], 



(IV -5) 



a 



K 



8.30 f . 2.2 -I 

r z 4n 1 ' n J 



(IV -6) 



The 4> dependence of Qr tan ^ reflects the use of both stereographic 
views in the determination of tank. 

2. Turbulence Errors 

The superheated condition required for track formation is 
obtained by rapid expansion of the liquid hydrogen prior to entry of 
each beam pulse. Because of turbulence following each expansion, 
tracks in the chamber are distorted, and curvatures obtained from 
photographs of these tracks can differ considerably from the curvatures 
of the actual ionizing particles . 

14 

The momentum -err or data of Cresti were taken as the basis 
for including turbulence effects within the total curvature error, 
Measurements of 386 beam tracks at 1.12 Bev/c momenta were made. 
The square of the standard deviation of curvature can be expressed as 



N n 2 



1 


V 


1 


lL\ 


N-l 


z_ 

i=l 


P. 


Ip) 




L 1 


\ / mean_ 



(IV -7) 



where (l/p) is the mean value of measured curvatures. Data were 

obtained by Cresti under the following conditions: (a) photographs were 

taken with 3 msec, light delay; (b) tracks were in the central region 
of the chamber (the exact form of turbulence is dependent on the 
position of the track with respect to the chamber expansion line); 
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(c) the average interval between measured points was 1.0 cm = Meas- 

urements were made on beam -track segments of varying lengths; three 

average track lengths were used to determine the R dependence of 
2 

0 ^ . The average lengths of the tracks considered were 6.1, 9.8, and 

19-8 cm. The three points of Fig. 5 were determined by use of Eq„ (IV-7) 

and the Cresti measurements. 

A fit of these points is possible with a quadratic expression 
2 2-2 

0 ^ ~A + Bx + Cx , where x = R and R, the half chord length of the 

projected arc. Such a fit is indicated by inspection of the points. The 

correspondence between the functional form of the internal error in 

curvature, G ^ in E^fJV- 6 ), and the last term of this expression made 

possible a natural division of the total error, i. e . , G ^ = (G„ + (a ' )^, 

ii t ix K Jx 

where G ^ and G ^ are the turbulence and internal-measurement errors 
, XV Jtv 

respectively. The dependence of C ! on N, the number of points meas- 
ured on a track, was taken from the internal-error expression as 
~ C ’ 15 

C(N) = (1 - ). For flat tracks with points taken at 1.0-cm intervals, 

one has N =2R, hence the expression 

o 2 = A + BR " 2 + 0.5 C' R " 4 [R " 1 ( 1 - 0.75 R -1 )] (IV- 8 ) 

XY 

was fitted to the points indicated. The coefficients A, B, and C ? 
required for this fit were 

A = 9.32 X 10“ 3 , B = 0.174, C' = 62.5. 



Repeated measurements of single events were also used to distinguish 

between internal errors and turbulence errors. 

In the kinematic analysis of events, Section V, the quantity 

£./a> , where £. = (x. - x™)/ a., ( 1 = 1 , 2 , ••• 9 ), is a measure of the 
i i ii l 

shift made in the ith variable to obtain a fit to a designated 
interpretation. The distributions of and £ 3/0 £ over a 

number of different events are a measure of the fidelity of the total - 
curvature error formulae when applied to the outgoing negative and 
positive track respectively e Histograms of these distributions are 
shown in Fig. 6 . 

The distributions of ^/q ^ over many different events 
complement the repeated-measurement technique, and these distrib- 
utions were used in the adjustment of the err or -formulae coefficients. 
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Fig. 5. Curvature error from study of 386 beam tracks 
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Fig. 6. normalized adjustment of variables, 84 elastic 

events ** (with stopping proton). 
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The error formulae as programmed for normal nonstopping tracks 
are 

a R = [AtB/R Z 4 C/R 4 ]^ fBev/c) ~ 1 (IV-9) 

where 

A = 18.64 X 10" 3 
B = 0.226, 

C = (62.5/N)(l - 1.5 N' 1 ), 

R = projected chord length (cm), 

N - number of points measured in View I; 

? 2 — 

O ^ - [CJ ? ^ -t ] 2 (degrees) (IV-10) 

where 

2 1 — 

= 1.85 R N " 2 [1 - 1 / N ] , 

a "<j> 2 = (°- 188R ) 2 (ArB/R 2 ); 

and 

0 v = ^Ar- (1 + 2 j sin <j> | ) " 1 . (IV -11) 

tanK RN* 



The O . term in Eq, (IV- 10) reflects the effect of turbulence 

T 1 1 

on azimuthal angle <j>„ It can be shown that this turbulence error, O 



and the turbulence error in curvature, 
expression 



K 7 



are related by the 






Cr'^rad. ) = (R/305)a" K , 

where 305 cm/(Bev/c) relates curvature and momentum in the magnetic 

field, and O ^ = (A + B/R ) 2 , as specified in Eq. (IV-9). 

For protons that stopped in the chamber, the momentum (and 

hence the curvature) was determined from the range -m omentum 

0 273 

formula, p = 0.141 f * , where i is track length (cm). This 

expression, based on the range -momentum curves for liquid hydrogen, 
defines the proton momentum much more accurately than is possible 
when curvature measurements are based on parabolic fittings. 
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The range -momentum expression can be approximated by 
p - 0.141 i 4 for determination of curvature error associated with 
measurement of a stopping proton. With half projected chord length, 
R, and dip-angle tangent, tan X, as independent variables, differen- 



tiation of this expression gives 



1 



p = 4K 



R 



(Bev/c), 



, + (sin X O .) 

R cosh tan X. 

where i = 2R/cos\. For "flat'' tracks this expression reduces to 
0 K ~ KCT r / 4R. A value of 0.055 cm was taken as the error in one-haif 
chord length, 0 R . This value was the standard deviation in range for 
444 p ranges measured by Clark and Diehl, using the same meas- 
uring engine. The approximation programmed was 



O 



K 



= 0.014 K/R. (Bev/c)" 1 



(IV -12) 



Since the error in K for stopping protons is essentially 

internal-measurement error, 36 remeasurements of event No. 495,621 

were made to test this approximation. The arithmetic mean of the 

36 0^/s f° r stopping proton obtained by Eq. (IV-12) was 0 072 (Bev/c) \ 

The K~ distribution for these remeasurements, shown in Fig, 7, has 
^ ~ 1 

a standard deviation of 0.128 (Bev/c) . The disagreement between 
the calculated error and the error associated with a typical event 
indicates the dependence of range -measurement accuracy on track 
characteristics. The measurements by Clark and Diehl were re- 
stricted to flat tracks with angles of inter section- -at the vertex-- 
between 30° and 150°, and to the less heavily ionizing muon and electron 
tracks--that is, conditions not normally encountered with actual re- 
coil protons. These considerations suggested that the expression for 
error in K be broadened to 

0 T , = 0.024 K/R. (Bev/c)" 1 (IV -13) 

.K 

The typical accuracies of the two methods of measuring 
curvatures (momenta) are illustrated by event No. 494,808: 

track length (i) = 3.034 cm. ; 

X = 58.59°; 

momentum by curve fitting, p^ - 121.3 Mev/c; 
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Fig. 7. Curvature of stopping proton; 36 repeated measurements 
of event No. 495,621. Kj = 9.68/Bev/c; cr ^ = 0.128/Bev/c. 
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momentum from range, = 192.2 Mev/c; 

momentum from least-squares adjustment, p^ = 188.5 Mev/c; 
momentum -measurement error from O^. Eq. (IV-13), G =5.0 Mev/c. 
Here and p^ differ from p^, the "best* 1 value of the momentum, by 
2% and 36%, respectively. The adjustment of p^ to obtain a ’’fit 11 was 
0.74 O . The inaccuracy of p^ reflects the short projected length of 
track (1.63 cm) available for the parabolic fitting. 

3. I ncident-Beam Momentum 

The nominal beam momentum of 0.927 Bev/c was determined 

, yz 

from a study of pion trajectories through the collimation system. 

The absolute value of the mean pion-beam momentum was determined 
from an analysis of 77 elastic events in which the proton stopped in the 
chamber. The nominal beam momentum, with a standard deviation of 
0 .278 Bev/c, was used in the kinematic fit program discussed in 
Section V. This large, unrealistic error allowed unrestricted shifting 
of the beam curvature, , in obtaining a fit to the intr epretation. In 
effect, no exact beam momentum was assumed. The adjustment 



parameter 



■ 



X 1 “ X 1 



M 



where and x^ are respectively "best fit" and nominal values of 
Kj, and ~ 0.278 Bev/c, indicates the adjustment of 

required for each event to fit the interpretation. The histogram 
for the 77 events is given in Fig. 8. The mean of -0.0414 ± .005 
corresponds to a mean beam momentum adjustment to 0.939±.006 Bev/c. 
The mean beam momentum, from a study of three associated production 
events of the type tt + p -► A + and tt + p obtained during 

/ 15 

this exposure, gives p^ = 0.938 ± .002 Bev/c. Abeam momentum 
p i = 0.939 ± .007 Bev/c was used. The indicated beam spread, 

± 0.007 Bev/c, is supported by a detailed study of associated pro- 
duction at 1.12 Bev/c. ^ The beam trajectories through the magnetic 
steering magnets and through the collimation system were in all 
essentials identical with the 0.927-Bev/c experimental setup. The 
study of stopping protons, on the other hand, indicated only that the 
beam spread was less than ± 0.035 Bev/c, a width that reflects the 
large intrinsic measurement error of the method. 






- , M v* JU1 •, , 

*••««*•*> (i ~*r •»* • •'« * » r- «»<i 

*< . • »•*•»» 

•• 





-23- 



Fig. 




Adjustment of beam curvature for least-squares fit, 
0.414 ± .005; a t = 0.1727. 



Mean = 
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V. KINEMATIC ANALYSIS 
A. Genera l 

In the past, it was possible for the physicists using cloud 
chambers to keep up with the data accumulation by "hand analysis " 
techniques. The bubble chamber, with its relatively high density, has 
greatly increased the number of interesting events observed, and 
therefore has stimulated the development of rapid analysis techniques. 
In this experiment an event was observed with approximately one out 
of every two Bevatron pulses, or 300 events were photographed each 
hour. The analysis of the events photographed during one 24-hour 
period of operation would require several man-years of "hand 
analysis. " Rapid precision analysis techniques are obviously required 
if the rate of analysis is to parallel the accumulation of data. 

All computations described in this report were carried out by 
use of the IBM type 650 digital computer. The procedure followed 
in the analysis of two-prong events is outlined in Fig. 9. The meas- 
urement and event-reconstruction phases of analysis have been dis- 
cussed in Section IV. 

B . Computation of Measurement Errors and Coplanarit y 

The origin of measurement errors and the derivation of 
requisite formulae are discussed in Section IV. The computer is 
programmed to calculate the errors associated with all angle and 
curvature measurements. A section of this program computes the 
included angle between any two tracks; individual track errors are 
propagated to these included angles. 

The coplanarity of an event may be measured by the value of 
R = (n^ X n^b where rip n^, a re the unit vectors in the 

direction of the incident ir and the two outgoing prongs, respectively, 
at the vertex. Perfect coplanarity is, of course, indicated by R = 0; 
however, because of measurement errors, events that are actually 
coplanar may have values of R ^ 0. A coplanarity function F and its 
associated error, CFp, were used in establishing the degree of co- 
planarity of each event. Here R and F are related by 

R = F cos k^ cos k^ cos k^, 
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Average time : 



Average time : 



Average time: 



Average time : 



Average time : 



mm 



2 min 



0.2 min 



3 5 min 



0.3 min 



Fig. 9. Flow diagram for analysis of two prong event; computer 
time indicated. 




Since tracks dipping or rising vertically in the chamber are not seen, 
cos \ is never zero, and R = 0 if and only if F = 0. The error, <7p, 
is determined directly from the measurement errors computed for 
each track. The quantity iFl/ctp indicates coplanarity within meas- 
urement errors. Figure 10 is a histogram of | F 1 / O' values for 
elastic and inelastic events. For the inelastic events a broad 
distribution of |F]/ is expected, since the lack of coplanarity is 
determined by the direction and momentum of the neutral particles, 

C . Graphical Analysis of Inelastic Events 
1 . Ellipse Plots 

In the following discussion, inelastic events 

- , , 0 
7T + p TT T p T TT 



and 



tt + p“*it + t r +n 
are represented symbolically as 

(Mj) + (M 4 ) - (M 2 ) + (M 3 ) + (M 5 ). 



Subscripts on kinematic variables indicate the particle to which the 
variable applies. Unstarred quantities and starred quantities pertain 
to the laboratory and center -of -mas s (c.m.) systems respectively. 

In the laboratory system that component of the momentum of the 
neutral particle, M^, that is perpendicular to the plane defined by the 
two charged particles introduces a deviation from coplanarity and 
prevents the use of a geometrically simple dynamical representation 
in this reference frame. In the c.m, system the three final-state 
particles are coplanar, and, in this frame, ellipse plots may be used 
in the kinematic analysis. The ellipse plots for analysis of the two 
inelastic reactions are shown in Fig. 11. 

The properties of each ellipse are established by the Lorentz 
transformation from the barsy system, (primed quantities) in which 
the virtual mass of particles and is at rest, to the c.m. system 
(see Fig . 1 Z ). 
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Fig. 10. Coplanarity of 250 elastic events (solid line) and 250 
inelastic events (dotted line). 
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Fig. 12. Ellipse parameters, with diagram 
from barsy to center -of -mass system. 



of transformation 
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Thus for particle No, 3, (p 1 ^ = -p ! ^ = p 1 ; 6 i = 0 =tt - 0^), one 
can write 

/ * *' 

'p ~ cos 0 




only the longitudinal component of the momentum -energy vector being 
transformed. The quantities y and q refer to the motion of the barsy 
relative to the c m. system. 

The origin of the ellipse is given by considering the special 
case in which the particle is going off at tt/Z in the barsy, 

# * _ 

p ^ cos 0 ^ - t) E y 

The major axis is the difference between the momentum vectors 
p ^ for the two special cases with 0 1 = 0 and 0* = tt, i. e, , Zy p^. The 
minor axis is given simply by Zp j ^„ 

Computation of the ellipse parameters necessary for construe- 
tion of the ellipse plots was programmed for the IBM 650 by the 
following procedure. 

(a) Calculate the virtual mass of particles No 3 and No. 5, 
>rt 35 = [E * + M 2 - Z E E 2 ]2 , 

* 

where E is the total center -of -mas s energy, 

* _ # 2 Z - 

(b) Calculate p ^(max) = [(E^ (max)) - J 2 , 

where 

E ^ (max) - — ^ [E 2 + M 2 2 - (M 3 + Mg) 2 ]. 

* 

(c) Compute the ellipse parameters for increments of p ? in 
the interval 0 < p^ < p^> (max), using the following expressions: 

p ' ~ Z>^ 35 1^35 + M 3 r M 5 ” f? ^35 ‘ M 3 * M 5 ) ^35 + M 3 M 5 ^ 

C^ 35 - M 3 + M 5 )]2, (Bev/c) 

f, E' 3 - 0.5 p z [l + (M 3 2 - M 5 2 )/ 7n 2 35 ] ) Bev/c 

V = (E* - E* 2 )/^ 35 ; 



where 
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In both inelastic cases considered, the negatively charged 
particle is a n , and the positive particle in the final state characterizes 
the event. 

Z. Transformation to Center -of -Mas s System 

a. Rotation of Chamber Coordinate System 

For the graphical analysis, a transformation of the measured 
momenta and angles of the prongs into the c.m. system is required. 
Propagation of the errors associated with these measured quantities 
into the new reference frame is also necessary. 

Error propagation into the c.m. system was simplified by 
first rotating the measurement coordinate frame so as to bring the 
incident beam direction into coincidence with a transformation axis. 

In the coordinate system of the chamber, System I, Fig. 13, the 
direction vector of an incident tt is nearly parallel to the positive y 
axis. System I was rotated to bring the beam direction into coincidence 
with the positive x ! axis of System II. 

The beam track in System I has polar and dip angles (<|>^»kj) 
defined, in terms of direction cosines, by 

Z — 

sin X., = n, ; cos X, = 4* [T - n. I 2 ; 

1 i z i i z 

coscb, = n, /cos L ; sincb, = n. /sink.. 

I lx 1 1 ly 1 

It can be shown that rigid rotation of the coordinate System I 
through an angle of <\> ^ about the x axis, followed by a rotation through 
an angle of -X^ about the new y axis brings about the desired align- 
ment of the incident track and the +x’ axis of System II. Thus, a 
prong n represented in System I by 

n = cos X cos <(> l + cos X sin <|> j + sin X R 
will, in System II, have the form 

n = sin a cos (3 i' + sin a smp j 1 + cos ak‘, 



with 



/sin a cos (3 
sin a sm(3 J = M(4>j,X^) 
y cos a 



cos X cos <jv 
cos X sin <j> 
sin X 
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Fig. 
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The real orthogonal 






' COS X^ cos 4> 1 



cos X^ sin^ 




cos 4> 1 




In System II the direction of the kth emergent prong is 
specified by 

cos a^_ = cos X^ sin X^_ - sin X^ cos X^_ cos (c[)^_ - c|) ^ 
sin a^_ = + \/ 1 - cos^ a^, 

cos p, = — r - [cos L cos X, cos(4>, - 4> 1 ) + sin X, sin X, ] , 

1C Sin CL 1 iC K. 1 1 K. 

k 

sin(3 k = slrTa, [cos X k sin (4, k - <(>!)]. 
k 

It can be shown that the errors in a and (3 associated with 
errors in <J) and X are, for the kth prong, 

\ 2 ■ A il t«v 2+! V 2! + A » ( V 2 + A 3k %> 2 ’ 



° P 2 - B u [«y 2 + <y 2 i + B zk <y 2 +B sk <v 2 - 

where the coefficients A^* are g^ ven b y 

A, - = sin 6 , sin X, , 
lk r k 1 

A 2 k = cos P k - 

A 3 k = cos X^ cos X^_ + sin X^ sin X^_ cos - <f> ^ ), 

B lk = ~T 2 ~ COsX k A 3 k ’ 

Sin “k 

B 2 k = isrs. cos “k Sin( 3 k s 

k 

„ _ sin Xj sin(3 k 

3k 



sin a k cos X k 
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b. Center -of -Mas s System Quan tities 

Following the axis rotation the incident -partic le track coincides 
with the + x 1 axis of System II, and this axis is the transformation axis 
and the reference axis for angle measurements; 0^ is the laboratory 
space angle for the kth prong. 

For the kth prong, we can write 




and the center -of-mass momenta necessary for graphical analysis 
are given by 

p k = p k [ sinZ \ + Y 2 ( cos - i- ) 2 ] 2 > (k = 2, 3) (v-l) 

where 

Pk 2 Z 

(3 1 = n™ ■— ~i< — , cos 0, - sin a, cos (3, , sin 6, = 1 ~ cos 0, . 

k n/ p2 +m 2 k k k k k 

k k 



The angular separation of the two emergent prongs in the c.m, 

* 2 

system, ® 2 3 * can ca ^ cu ^ a ^ ec i from the invariance of ^ ^ , i. e. , 



^23 = M 2 2 + M 3 + 2 E 2 E 3- 3 



P2 P 3 cos ®2 3 



hence 



2 2 t 1 * ^ 

= + M^ + 2 E^ E^ - 2 p ^ p ^ cos 6 > 

a * -1 ( E Z E 3 " E 2 E 3 + P Z p 3 COS 0 Z3 

0^3= cos 



5§C * 

p 2 p 3 



with cos = sin a ^ sin cos (p^ - p^) + cos cos 



(V-Z) 



and 



* 



E k = / 



*2 

Pk 



+ M 



From Eq. (V-l) it can be shown that the standard deviation of 
the error associated with the kth prong's momentum (c.m.) is 
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where 



Z - 



- -* [< c kV +(D kV» 



k Pi 



„ , p k ] »» M k - 0 

k ~k E kP k v k ' p k 



D k = P k [sin 6 k cos 0 k (l-7 Z ) + ^ sin , 



and 



|OT I (sina k sin Pk |Z ^ + (cosa k co S p k ) 2 a 2 P 



* 

The error in 6 a function of the six variables 
(p , p^> a 2 > $ 2 * a 3 , P 3 ) that define tracks No. Z and No. 3, was 
obtained by evaluating the required partial derivatives of Eq. (V-Z) 
as difference quotients rather than programming in their analytical 
expressions . 

The center -of -mass quantities are, of course, dependent on 

the reaction assumed, i. e. , the mass assigned to particle No. 3. The 

computer program first computes the quantities using = M ; then 

P 

the calculation is repeated with = M^. 

For each event the two inelastic reactions are considered in 

* * 

turn; and the areas, defined by p^ , cos 0 ^, an< ^ their associated 

errors, are plotted on the appropriate set of curves. For one of the 

% 

two reactions the terminal point of p~ will, ideally, lie (within errors) 

6 * 

on the ellipse specified by the value of p^ , thereby characterizing the 
reaction. A typical event has been plotted in Fig. 11. The values of 

p->, p~, and the standard deviation O * , O * are used in determining 

23 P 2 P 3 

the best graphical fit. The quantity to be minimized in each instance 




= p*, O. = O * ( i = Z, 3). 

*11 p. ' 

*1 



where 



x. 

1 
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One can write 



Z 

* (min) 



(x 2 - x 3 )' 

tt 



The values of \ (rnin) ^ or t ^ Le ^ ne ^- astlc interpretations were compared 

for determination of the correct graphical interpretation. For the 

event plotted in Fig. 11, the missing neutral is obviously a neutron; 

Z Z 

X tt-o/x n “ 99. A detailed discussion of the graphical separation of 
inelastic events is not given. In summary, separation by this method 
was identical to the separation obtained by the computer program 
discussed in the following section. 

D. Computer Program for Kinematic Analysis 
of Two-Prong Events 



1. General 

Kinematic analysis of two -prong events — i. e. , events having 
one incident tt track and two charged outgoing tracks --was programmed 
for the IBM-650. The input data consist of the measured values of 
nine dynamical variables x^ (specified below) characterizing each 
event, and for each measured variable its associated standard 
deviation, Ch . Output data consist of (a) the least-squares adjusted 
value of each of the nine variables required to make a "best fit" to 

each of three possible interpretations, (e. g. , elastic, and inelastic 

0 , Z 

with missing tt or neutron); (b) the value of x a nd a number of 

allied statistical consistency tests; (c) the momentum and trajectory 

of the missing neutral particle (rr^ or neutron). 

Z. Choice of Dynamical Variables 

The nine dynamical variables x^(i = 1, Z, ••• , 9) consist of 

three curvature variables (j = 1,Z,3), 



K. = -*! 

J I p • I c o s X. ■; 
J J 



(Bev/c) 1 



where 



± = 



+ 1 for incident track, track No. 1 

-1 for outgoing prongs, tracks Nos. Z and 3 










Wj* 4 * ' - * 
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the three dip tangents 

tan Xj, ( j = 1,2,3); 

and the three azimuthal angles 

4 >j (j - 1,2,3). 

These variables were chosen because of correspondence with 
quantities most directly measured from the stereo projections on a 
scanning table. Curvature was chosen in preference to the more 
familiar momentum variable to eliminate the following problem: the 
use of momentum as one of the fitting variables introduces an artificial 
mathematical discontinuity at high particle momentum. For example, 
two beam tracks of momentum 4.0 Bev/c and opposite charge are 
indistinguishable by present techniques. Their momenta, however, 
are two large numbers of opposite sign, whereas their curvatures are 
essentially zero in both cases. At lower momenta this difficulty again 
shows up as a skewness in momentum error. 

3 . Interpretation Define d 

The general two-prong interaction may be written symbolically 
as 

(Mj) +(M 4 )- (M 2 ) + <M 3 ) + (M 5 ), 

where 

= incident particle (track No. 1) 

= struck particle, 

= track No. Z particle, 

= track No. 3 particle, 

M<- = missing neutral (not used if elastic). 

An interpretation consists, by definition, of two statements: 

1. Whether the event is to be treated as elastic or inelastic; 

Z. The set of masses M^, M 3 , -- and, for inelastic inter- 



pretations, 
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4 * Kinematic Equations of Constraint 

For any given interpretation the momentum and energy (-)F^ 
{\= 1 , 2 , 3, 4) of the missing neutral, if any, can be written as 

3 

F = (-) Y. cos (V-3) 

j=l J J 

3 

F z = (-) ^ sin<t>./K, (V -4) 

3 

F 3 =(-) Y. tanK/K, (V-5) 

j=l J J 



F 4 = (-) (M 4 + E l - E 2 -E 3 ) a 



(V-6) 



where the energy of the jth particle is computed from 

(1 + tan^\.) , 



E. = 
J 



J 



K. 

J 



+ M 

J 



with M the mass of the j_th particle. Here M^fBev) is the mass of 
the struck particle. 

The kinematic constraints can be summarized into two groups 
according to whether the event is treated as elastic or inelastic: 

1 . Elastic: Impose four contraints, F^(x) = 0 for \ = 1 , 2, 3, 4. 

2. Inelastic: Impose the single constraint, £(x) = 0, 
where t, is defined by the expression 

(M, - ]Q|^ Q £0 

i= 5 

[m 5 + IQ I 2 Q < 0 

where 3 

Q = (-) F v F x - F/ - I F 2 . 

1= l 



_1 

For Q positive, lQ| 2 is the apparent mass of the missing neutral. 

The quadratic dependence of £, on F^ makes it possible to 
reach a least -squares fit for which the energy of the missing neutral 

is negative, i.e., F^>0. These erroneous solutions are fortunately 
infrequent and easily detected. 
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The particular kinematic constraints selected above are some- 
what arbitrary. For example, one of the elastic constraints could 
have imposed the coplanarity restriction. Almost any function 
G(M^-Q) such that G(0) was zero could have been used for the inelastic 
constraint. The kinematic constraints cited above were selected first 
for their ease of computation, and second for their essentially linear 
dependence on the dynamical variables in the range of physical 
inter e st. 

5. Definition of a Kinematic Fit 

In general the measured values of the nine variables of 

an event do not satisfy the kinematic constraints. Therefore a least - 
squares adjustment of this input data is required with the measured 
quantities, x^ s going to the adjusted values, x^: x.^ -► x^ (i = 1 , 2 * • ■ 9) . 

A least-squares fit to a given interpretation is specified by a 
set of dynamical variables x^ satisfying two conditions: 

(a) The constraints are individually equal to zero, i. e. , 

F^(x^) = 0, for all constraints X = 1, c. For inelastics 
F^ is set equal to £> and c = 1. 

(b) The chi-square quantity, 




is at a minimum (or stationary) point for all variations, 

6 x. , about x. consistent with the kinematic constraints. 
i i 

An alternative, but equivalent, formulation using Lagrangian 

multipliers proves more mathematically tractable. This choice 

replaces (a) and (b) by this requirement: 

The modified chi-square quantity 

c 



M(x., a^) = x' 



+ 2 



I 

X=1 



w 



(V-7) 



is minimum (or stationary) under arbitrary variations, 6 a^, 6x^ of 

all variables about their stationary values, a^, x_ 













1 . 
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6 * Statistical Estimates of the Quality of the Fit 

Achievement of a 1 ’fit 11 is a purely mathematical achievement 
in programming; in no way does it assure that the selected interpretation 
is valid. Therefore, a statistical estimate of the quality of the fit 
to a given interpretation is required. 

The adjustment parameters are normalized by introducing the 
quantities 

g. = (x. - X. M ) /O. (i =!,-•• 9). 



The foremost quantity for the statistical estimate is, of course, 



Z 

X = 




evaluated at the "fit." If y ^ is in the doubtful zone, being neither 
unduly large nor sufficiently small, other tests can be quite useful. 

In particular a distinction is made between: (a) "errors' 1 
from inherent inaccuracy of measurement and from turbulence, as 
discussed in Section IV; and (b) "bloopers" caused by miscoding, 
measurement of wrong fiducials, or jumping to a neighboring track, 
or digitizer errors. In addition, it is helpful to distinguish three 
statistical estimates that were used. 

(i) Globa l: the tests < 3. These tests are dependent only upon 

x at fit, the measured x.^, and the standard deviation in the variable, 

l l 

a . 



(ii) Loca l: the tests | O ^ | < 3. These tests are peculiar to the 
interpretation fitted. Here 1 (7^ is the standard deviation in £. itself. 
Such tests exploit the derivatives of the constraints at fit and are often 
more sensitive than the global test, provided the interpretation is 
correct. 

(iii) Systematic: correlation of curvature and azimuthal angle errors, 

incorrect nominal beam momentum, and neglect of dE/dx corrections 
in curvature measurement. Corrections of this nature were made by 
hand calculation where required for resolution of otherwise ambiguous 
interpretations . 



. 
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? • Trajectory of the Missing Neutral Partic le 

The variables appropriate to the missing neutral track-~i.e. » 
curvature, ; dip tangent, tan ; and azimuthal angle, 4 '^ '-are 
calculated for each inelastic interpretation. When a fit has been 
obtained, these quantities are specified by 



K 5 = (p 5 cos \ 5 f 1 = (F } 2 4 F*) "2, 
tan X- 5 = F 3 /fF 1 2 + F^)* = FjKg, 



(Bev/c)' 1 , (V-8) 

(V-9) 
(V-10) 



<J)j- - arc tan (F^/F^), 

where ~F p -F -F^ is the momentum vector of the missing neutral, 

Eqs . (V-3), (V-4), and (V-5). 

8, Mathematics of the Fitting Procedure 
a. Linear Constraints 

The procedure introduced here was developed by Frank T. Solmitz 

16 17 

and Horace Taft. * Let us examine the fitting problem for the 
special case, in which the constraints constitute linear relations 



between the x^, 



F \< x i> = *X + S- F \i x i’ 
1= 1 



(v-li 



where F^. = 3F^/9x., We now require the (9+ c) partial derivatives 
of M, Eq. ( V — 7 ) to be zero: 



3 M » e V TT 

_ — = 0: § = - 2_ a. F. . a . 

ox. i q- X Xi i 



(V - 1 2) 



|^- = 0: F. (x.) = 0. 

8 a. X i 



(V- 1 3) 



Substituting Eq. (V~12) into (V-13), we get c equations for the a^'s: 



La 



H,„ = F, <x M ), 



X Xp X 



(V-14) 



where is the symmetric matrix, 



H 



Xp 



Y (F, . cr )(F a ) 
2 -, ' Xi i ^ pi i 



i=l 



(V-15) 
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Hence we have 



a. - 2 H 



-1 



F.., 



\(J. (X 



(V- 16) 



which, when substituted back into Eq. (V 12), gives the adjusted values 

of X. . 

1 

b„ Nonlinear Constraints 

Solmitz has suggested the adaptation of this method to the 

solution of the fitting problem for nonlinear constraints.^ For a point 

v 

in the neighborhood of the fitted values x^ the nonlinear constraints 
can be approximated by the linear expansion 

F x ( Xi ) = F x (x. V ) + aF/ax i l x v (x. -x. V ). (V-17) 

i 

The method of successive approximation, using equations paralleling 

Section VD8(a) is then used to derive better estimates, x^ V , of the fit 

variables x. . 

i 

The fundamental Eqs. (V-12) and (V-13) carry over directly; 
the quantities F^ must be here reinterpreted as partial derivatives: 

F u = hr F X ■ < v - 18 > 



Given the vth approximation (x. v , a . v ) to the solution of the fundamental 
Eq . (V - 1 2) and (V-13), one deduces the next iteration (v + 1 ) fr om a 
first-order expansion in (x^ -x^ v ): 



f. V+l V V+l ^ V 

^ = - L- a. 



F (7., 

Kl 1 



(V-l 9) 



F v v + i (C +1 -C)Fj <y ■ = 0 . 



i=l 



ki i 



(V-20) 



v+l 



Introducing the two quantities generalizing Expressions 

(V-15) and (V-16), and the modified constraint function f^ , one has 

H V = (F» V a.)(F V a.), (V-21) 

k(x < ~r ki l pi i" 

V+l V /tjV.- 1 v 

+ = L (H kix V ’ 



M- 



(V-22) 
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, V „ V . V > V V 

f \ = F K + 2_*i F Xi a i : 



(V -2 3) 



t 



v+1 



v+1 



is obtainable by substitution of into Eq. (V-19). 

This iteration process can be initiated with "zer oth-appr oxi - 

n 1 (0) M . t (0) n 

mation values x = x. , i. e. , = 0. 

l l i 

9. Statistical Correlation of Adjusted and Measured Variables 

M 

The measured variables x are assumed to be uncorrelated: 

i 



* M * M c 2 
Ox Ox . = 0 . . o' . , 

i 3 iJ i 



where a. 



M 



is the standard deviation in the measurement of x^ . However, the 
fitted variables x^ are correlated as a result of the fitting process. 
The normalized displacements £ (i - 1, ••• 9) are also correlated. 

Evaluation of the fit achieved for any interpretation is possible 
only if we establish two references: (a) the standard deviation , 

cr p , of itself; and (b) the average value of y to be expected with 

i 

c constraints, using the fitting procedure. 

2 

These results can be deduced as follows. At \ minimum the 
normalized displacement ^ is given by 



t = - 



i 



f,, H 

(J- (J.V 



-1 



(F . (T ), 

' VI 1 



and accordingly the variation 6 is 



st 



-I 



[±,V 



6f L1 H ” (F . (J ). 
M* [xv 1 vi i 



16 



Using the result of Solmitz , valid only to the extent that the region 
Z 

near \ minimum can be treated as linear, 



5 f, Sf = H , 

M* v [iv 



a 



t 



= (5 U‘ 



= I 



v=l 



(F . a.) H (F . a.). 

[XI 1 [XV ' VI 1 



then one deduces 

































4 |< M 
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Not all standard deviations a * are equal to unity, because 
in a least-squares fitting those dynamical variables least well known 
are, in general, adjusted (shifted) by a greater number of standard 
deviations in the fitting process. The digital computer was pro- 
grammed to evaluate and print out the adjustment variables ^ 
which have in principle a unit-width normal distribution. 

The mean value of can now be established as a special 

result : 





c 

. i 

1, [X, v=l 



(F 



. ct.)(F cr . ) H 
pi i vi r 



-1 

pv 



= c, the number of constraints. 

In summary, x^ should average to 1 and 4 for the inelastic 

2 

and elastic events, respectively. The x distributions for the inelastic 
and elastic events are shown in Fig. 14. 

1 0 . Comparison of Measured and Computed Quantitie s 

Figure 15 shows a photograph of a ir + p tt +tr^ + n 
interaction. A schematic diagram of the event is shown in Fig. 16, 

This is a typical two-prong event if we consider only tracks Nos. 1, 2, 
and 3. (The kinematic program considers only these tracks.) In this 
event, the missing neutron (neutral track No. 4) strikes a proton, 
track No. 5, and the neutron momentum and trajectory can be calculated 
from the direction and range of this stopping proton. This fortuitous 
circumstance permits the following tests of the kinematic program, 

(a) Does the program, utilizing only measurements of tracks Nos. 1, 

2, and 3, achieve a n fit M for the neutron interpretation? (b) Do the 
computed trajectory and momentum of the neutron agree with the 
measured values (based on the vertex, the first bubble of the proton 
track, and the range of the stopping proton)? Table III summarizes 
the comparison between measured and computed quantities. The 
large measurement error associated with the dip angle of the neutron 



I W«| , i 

-< • ’ . t • r N«i 







NUMBER OF EVENTS NUMBER OF EVENTS 
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Fig. 14. (Above) 
events 

(Below) 



X^ inelastic interpretation, 250 inelastic 

2 , 

X elastic interpretation, 250 elastic events. 



-46- 




ZN -1941 



Fig. 15. Photograph of tt + p -► tt + tt +n event (neutron 
scattering) in the bubble chamber. 
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MU-15,340 



Schematic diagram of tt + p -► tt 



+ TT 



+ n event. 



Fig. 16 
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is caused by the shortness of the track and the fact that only two points 

were available for measurement . The measurement errors for the 
-f 

it (a typical track) are included in Table III. 



Table III 



Measured and adjusted quantities 






Value from 


Measurement 


Least-squares 


Partic le 


Quantity 


measurements 


error 


value 


proton, 


range 


7.178 cm 


— 




No. 5 


P 5 


0 .243 Bev/c 


0.002 Bev/c 




neutron, 


P 4 


0.256 Bev/ c 


0 .036 Bev/ c 


0.267 Bev / 1 


No. 4 


+4 


95.1° 


4.2° 


93.1° 






19.3° 


38.0° 


25.4° 




length 
(meas . ) 


0 .325 cm 


— 


— 


— 


x el 






>99.9 




* 2 „ 






2.0 


TT - *" 


P? 


0.388 Bev/ c 


0.025 Bev/c 


0.372 Bev/ 


No. 3 


*3 


74.7° 


0.2° 


74.7° 






11.9° 


0.5° 


12.0° 




length 
(meas .) 


13.138 cm. 


— 


— 
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VI. CLASSIFICATION OF EVENTS 
All events were first examined for the value of the coplanarity 
function F/cTp, and catalogued as coplanar (F/a < 4) or non- 
coplanar (F / a ^ >4). 

The noncoplanar events were considered to be definitely in- 
elastic, and were subdivided as follows: 

Z Z Z 

(A.l) Strongly separable by x. * i.e. , x X n ^5, or 

vice versa. Final classification of these events corresponds 
2 

to the smaller x • 

Z Z 

(A.Z) Tentatively ambiguous, i. e. , both x ^0 an< ^ X n ^ 5. These 
events were examined for delta rays and ionization of the 
positive track, and kinematically examined on the Pq, p^ 
plot described below. 

(A. 3) All others - not having satisfactory x^ for either the tt® or 
neutron interpretation. These events were examined on the 
scanning desk for measurement and (or) coding errors, 
remeasured or adjusted by template fittings, and finally 
forced into category A.l or A.Z. 

The coplanar events, | F / a ^ I < 4, were divided into classes 
in an analogous way: 

2 

(B.l) Strongly separable as elastic, by x 

(B.Z) Elastic or tt (ambiguous); neither ionization nor delta rays 
can possibly resolve the ambiguity. 

(B.3) Elastic or neutron (ambiguous); possible to resolve by 

identification of positive track, on the basis of ionization and 
delta rays. 

(B.4) Neutron or tt® (ambiguous); separable as class A.Z. 

(B .5 ) No fit achieved; events remeasured. 

All the events with |F/c7p,|< 4 and fitting only inelastic 
interpretations were re-examined for gross measurement errors. 

These errors can induce a spurious or neutron fit. 

Positive -particle identification, based on ionization, was used 
in the classification of otherwise ambiguous events. The incident tt 
was used as the local minimum -ionization standard. (Calibration of the 
chamber for bubble density as a function of ionization was not feasible 
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for two reasons: (a) chamber sensitivity was kept high to detect neutral- 
particie end points; (b) superposition of bubbles - -especially for steeply 
dipping tracks --makes bubble counting unreliable. ) Low-momentum 
protons were readily identified in the chamber. 

It is quite possible for an inelastic event to be kinematically 
ambiguous, even in the limit of exact momentum and angle measurements. 
To see this let p^_ be the observed momentum of the positive prong, and 
let Pq be the calculated momentum of the missing neutral (tt^ or neutron). 
Normally the energy available to the positive and neutral particles, as 
calculated from the conservation of energy and the known masses of 
the incident and struck particles, allows assignment of only one mass 
to the positive particle and one mass to the neutral particle (either 
proton mass and t r mass, or tt^ mass and neutron mass). However, 
to the approximation that the small n-p mass differences can 

be neglected, kinematic separation fails for p^ = p + . Although such 
events are kinematically ambiguous, they may, of course, be classified 
on the basis of delta rays or ionization along the positive track. 

Figure 17 shows the Pq, p + distribution of a group of 26 events 
of class A 2. The ambiguous p^ = p + line is broadened by the typical 

measurement error of p ; the error increases with increasing positive 

+ 2 
momentum. For the most part, those events with good \ separation 

have the better interpretation (tt^ or n) outside the ambiguity region. 



P o (Mev/c) 
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P + ( Mev/c ) 

MU-15, 341 



Fig. 17. Scatter diagram of positive and neutral particle 
momenta (26 events). Kinematically ambiguous region 
delineated. 
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VII. RESULTS 
A. Cross Sections 

The total number of interactions, corrected for over -ail 
scanning efficiency and detection dependence upon the orientation of 
the plane of emergent prongs is 1430 ± 43. After correcting for the 
small-angle elastic scattering, the total number of events occurring 
in the acceptance volume of the chamber is increased to 1452 ± 48. 

The total number of interactions and the integrated pion path 
length in the acceptance volume provides a direct measurement of 
the total cross section. Beam tracks crossing the acceptance volume, 
as projected in one view, were counted in every tenth frame of film. 
Based on this beam count and the average track length across the 
acceptance volume, a total pion path length of 7.76X10 cm± 1.1% 
was computed. This length includes corrections for muon contam- 
ination (1.0 ± 1.0%) and electron contamination (1.8 ± 1.2%). The 
7 r -p total cross section found is 53.3 ± 2.4 millibarns. The cross 
sections for ail zero-, two-, four-, and six-prong and V-particle 
events observed are listed in Table IV, together with their relative 
frequencies . 

Table V gives the relative frequencies and cross sections for 
the elastic-scattering and single ~pion-pr oduction reactions. 

B . Elastic Differential Cross Section 

Figure 18 shows the elastic 7r~-p differential cross section. 

This differential cross section displays the narrow diffraction peak 

in the forward direction, and after dropping to a distinct mimimum 

from 66° to 108° (center -of -mas s system) shows a pronounced hump 

in the backward hemisphere. The backward hemisphere contains 

35.6% of the events. The differential cross section at 0°, do (0)/6f2, 

extrapolated from the curve is 15 ± 4 mb/sterad. Cool, Piccioni, 

1 9 

and Clark give a value of 9.4 (center-of-mass system) on the basis 
of a 46-mb total cross section. Adjustment of the imaginary part of 
the forward-scattering amplitude for the larger total cross section 
of 53.3 mb yields a prediction of 12.6 mb for the elastic differential 
cross section at 0°. 
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Table IV 




Phenomenolog 


ical classification 








F raction of 


Cross section 


R eaction 


No, of events 


total events 


(mb) 


Zero prong 
Two -prong 


255 ± 9 
1182 ± 36 ’ 


0.176 ± 0.006 
0.814 ± 0.024 


9.4 ± 0.3 
43.4 ± 1.3 


F our ~pr ong 


12 ± 2 


0.008 ± 0.001 


0.4 ±0.001 


Six -prong 


0 


0 


0 


Associated 


3 ± 1 


0.002 ± 0.001 


0.1 ± 0.001 


production 








Includes ZZ ± 5 


❖ 

events in the interval 1 ^cos 9 

7 r 


>0.95 elastic 


differential cross section, Fig. 


18. 




Table V 




Elastic and sin 


tgle-pion production 






No. of events 


F raction of 


Cross section 


Reaction 


analyzed 


analyzed events 


(mb) 


Elastic 


261 


0.483 


21.0 ± 1.1 


“ . . 0 
t r 4* p 4* 7T 


80 


0.148 


6.4 


= , + , 

7 r + 7T + n 


199 


0.369 


16.0 


All other 
interactions 


— 


— 


9.9 


T otal 


cross section = 


53.3 ± 2.4 
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Fig. 18. Elastic 7T-p differential cross section at 0.939 Be v/c. 
Events corresponding to scattering of less than 18° in the 
center-of-mass system are omitted. 
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C. Distribution of Charged Particles in Single -Pion Production 

The momentum distributions of the tt and proton from 
tt -t-p-^TT -tp-f-TT^ reactions are given in Fig. 19. The momentum 
distributions are compared with those predicted by phase -space 
considerations. The average momentum of the observed tt is some- 
what higher than the average value predicted by phase-space arguments, 
and are consistent with the momentum distribution anticipated from the 
isobar model. ^ 

In Fig, 20, the momentum histograms of the charged pions 

— -f- 

from the reaction tt + p tt t it 4 n are drawn. Again the momentum 
distributions are compared with the phase-space distributions. A 
marked departure well outside statistics is seen; the tt^ peaks at a 
momentum of 125 Mev/c greater than center momentum of the tt 
distribution. 

No attempt was made to identify multiple -pion -production events 
during this experiment. The number of four -prong events observed 
(0 8 ± 0.1% of total events) indicates that pion multiplicities were 
infrequent at this value of incident tt momentum. The smgle-pion- 
production events analyzed in this experiment provide a reference 
point for analysis of single -and multiple -pion production at higher 
values of momenta, specifically the 1.12 and 1.24 Bev/c data now in 
the initial stages of analysis. 













. 









w r**m? II* * 





MU-15,343 



. 19. The momentum distributions of the protons and tt“ 
from the it" + p -» it" + p -{- tt ( - ) reaction. Expected phase 
space dis tr ibuti ons are shown for comparison. 
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7 T “MOMENTUM ( Mev/c) 

MU-15,344 



Fig. 20. The momentum distribution of tt's from the 

Tr~ + p -* TT“ + 7T + + n reaction. Phase-space distributions 
are shown. 
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